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Abstract.- We present a novel canonical description of the incompressible fluid dynamics. This de- 
scription uses the dynamical constraints, in our case reflecting incompressibility assumption, and leads 
to replacement of usual hydrodynamical Poisson brackets for density and velocity fields with Dirac 
brackets. The resulting equations are then known nonlinear, and nonlocal in space, equations for 
incompressible fluid velocity. 
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1 Introduction 

Canonical description of classical, compressible, isothermal fluid has been developed in the past 
[TJ [2l [3] for various purposes, for example description of superfluid 4 He 1 or in kinetics of the first 
order phase transformations [H |3] ■ An attempt to extend this formulation for the adiabatic flows has 
been proposed [3j and used to analyze dynamical properties of thermally driven flows. The isothermal 
flow canonical description can be generalize for the case of viscous fluids [5] within the framework of the 
metriplectic dynamics [5] . The Madelung representation for the wave function results in hydrodynamic 
like picture of quantum mechanics, where the "only" differences from Euler equations are hidden in 
the quantum pressure term, which is proportional to ti 2 , and in the quantization of the circulation, 
r = n(h/m). Apart of that the canonical description of the quantum fluid is then identical to that 
of the classical one. The dissipative generalization of the Schrodinger equation [7] also allows for 
metriplectic interpretation, which differs from the classical one [5]. 

The fundamental point in all of the above listed formulations of fluid dynamics is that the fluid 
density p is one out of the pair of canonically conjugated variables. For potential flows the other 
canonical variable is the velocity potential <f>. In case of general flow the two additional Clebsh 
potentials A, [i [3] in the velocity field representation v = —V<fi — Wp are canonically conjugated to 
each other. None of these descriptions can be applied to the case of incompressible fluids. 

The canonical description of incompressible flow is of considerable importance, for example in 
formulation of statistical mechanics of turbulent flow [5] . Many attempts to provide such a canonical 
formalism [10] failed to do so. The other important point is that real turbulent flow are hardly 
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incompressible, and the issues like compressibility corrections to scaling laws for turbulent flows are 
still open [TT] . 

In this paper we propose a novel formulation of a canonical description of the incompressible fluid 
based on the concept of the Dirac brackets [12]. The mathematical introduction to this formalism 
valid for general dynamic system subject to some set of constraints {9 Q = 0, a = 1 . . . N} can be 
found in |13) . Dirac bracket approach to description of incompressible membranes, within Lagrangian 
coordinates formulation of continuous mechanics of membranes, was given in [14] , We are unaware of 
any other application of that formalism to continuum mechanics problems. In separate publication 
we shall present other application of the Dirac constraints formalism in classical mechanics [15] . 

The Dirac brackets, for incompressible fluid, are presented, in what follows, within the Poisson 
bracket formulation of fluid mechanics [21 121 [5] , which avoids cumbersome introduction of the Clebsh 
potentials. Thus the state of fluid is described by specifying its density and velocity fields. 



2 Compressible fluid 

Consider infinite 3-dimensional volume of the isothermal fluid with density g(r, t) and velocity 
v(r, t). The Hamiltonian for such a system is given as: 

H{g,v} = f [pv 2 /2 + f(g)] d*r , (1) 



where f(g) is the fluid Hclmholtz free energy per unit volume, related to the fluid pressure by 

P=Qj- g -f{Q). (2) 

The Poisson bracket relations between fields g{r, t) and v(r,i) are [5]: 
{g(x,t),g(y,t)} = 0, 
{e(x,t)y(y,t)} = -Aj (x _ y)j 

{« i (x,t)y(y,t)} - i(x-y)-i-^(Vxy) t (x,(). (3) 

The continuity equation is obtained by evaluating the Poisson bracket {g, T~L}, and the Euler equation 
by{v,H} 

d t g{r,t) = {g,H} = -V-gw 1 

9 t v(r,i) = {v,H} = -v-Vv-(l/e)Vp( e ). (4) 

The above formulation of fluid mechanics can be derived from the least action principle provided we 
choose the proper lagrangian. As shown by Thellung [1] this lagrangian density is the local pressure. 

The incompressible fluid, although an obvious simplification, is adequate for all the flows when the 
local Mach number is small. The incompressibility condition then is that the density g(r, t) — go = 
what also implies that V • v = 0. Within the canonical formulation framework both these conditions 
are regarded as Dirac constraints 12J a (r, v,t) = 0, a = 1,2. Next section contains a brief overview 
of the Dirac brackets theory. 



3 Dirac Brackets 



The definition of the Dirac brackets we shall use in the following is a natural generalization for 
the original construction proposed by Dirac [12] and discussed in detail in [13]. When the physical 
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system with phase space V is subject to a set of constraints {6 a = 0} then its motion proceeds on 
a submanifold Pd5 = U^il 26 ^!®^-?) = 0}. If the Poisson bracket for two arbitrary (sufficiently 
smooth etc.,) phase space functions F and G was {F, G}, then the Dirac bracket C F, G □ is defined 
as: 



/V 



Ef,G3= {F, G} = - ^{F, e a }M ab {6 6 , G} 



(5) 



where M a (, is the inverse of the constraints Poisson bracket matrix W a b = {©a, Ofc}- 

The generalization of the Dirac bracket to the case of continuous variables, like in hydrodynamics, 
is straightforward. The sum over the indices a is replaced by sum and integration over the space 
variables and the inverse of the matrix W a b{r, r') = {9 a (r), 6 (r')} is defined as: 



]T / dr'W ac (r,r')M cb (r',r") = 5 ab 5(r - r") 



(6) 



Dirac brackets, given by Eq.© replace the original Poisson brackets in the equation of motion for 
the constrained system. Thus for a phase space function F the time evolution on the submanifold S 
is governed by: 



dF 
~dt 



(7) 



where T~L is system Hamiltonian. Next section will contain application of the Dirac brackets to the 
description of the incompressible fluid. 

4 Dirac Brackets for Incompressible Fluid 

The constraints used in constructing the incompressible fluid dynamics are: 



9i = q(t) - Qo = 
6 2 = V • v(r) = . 



The constraints Poisson bracket matrix W a b(r,r') can be evaluated using Eq.([3]), and it reads: 

e{r) c ' 



W o6 (r,r') = V;V J r 





s ij , 



i £ «^( V x v(r)) fc 



J(r-r') 



(8) 



(9) 



In the Dirac formalism one needs the inverse of the matrix W a b{v,v') defined in®. The matrix 
elements M a b(r, r') obey the set of partial differential equations, written explicitly in the Appendix A. 
Solving these equations we find matrix M a t(r, r') in the form: 



M o6 (r,r') 



M{G}, -G(r-r') 
G(r-r'), 



(10) 



where G(r — r') =| r — r' | /An is the Green function for the Laplace operator in the infinite volume, 
and 



M{G) = - J dx'G(r - x')V, 



1 



-V x -G(x'-r') x (Vxv(x')) 



(11) 
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5 Dirac equations of motion for incompressible fluid 

Using the definition and the explicit form of the Dirac brackets given in previous section and in 
the Appendix A, we first calculate the Dirac bracket C H □ . From the definition ([5]) , and ([7]) we 
obtain: 



dg(r,t) 
dt 



=E Q(t,t),n □= {g(r,t),n} -J2[ dzdz'{ e (r,t),e a (z)}M ab (z,z'){Q b (z'),H} 



(12) 



Explicit evaluation of the right hand side of Eq. (fT2"|) is a bit tedious, but using results from the 
Appendix A and (|23p one finds that it vanishes. Thus the continuity equation for incompressible 
fluid, within the Dirac formalism reads C 0(r,i),% □= 0. 

The algebra needed to derive equation of motion for the velocity field v is slightly more complex 
than these leading to the continuity equation. 

Following Dirac procedure we obtain: 



dv(r,t) 
dt 



=C v(r, t), H □= {v(r, t),U} - ^ J dzdz'{v(r, t), 9 a (z)}Af ab (z, z'){G 6 (z'), H] 

ab 



Evaluation of the right hand side of (JT3J), with use of expressions 
0v(r,t) 



gives: : 



Of 



v(r,t) x (V x v(r,i)) - V r 



J dzG(r - z)V z • {v(z) x (V x v(r, t))} 



(13) 



(14) 



Thus we have obtained nonlinear, nonlocal equation for the velocity field known from previous work 

GBH|. 

The above exercise in the Dirac brackets calculation provides a novel formulation of the Euler 
incompressible fluid. The viscous fluid equations can now easily be derived by replacing the Dirac 
brackets by the metriplectic brackets discussed in S\. We can also use the Dirac brackets as starting 
point in the perturbation theory in which compressibility corrections are calculated. To do so one 
formally associates small parameter k to the matrix elements M ao and expresses the Poisson brackets 
by the Dirac one. To the first order in k the expression is identical to that in ([5]) with reversed role 
of the Poisson and Dirac brackets. 

In conclusion we have shown in the above that the Poisson brackets formulation of the fluid 
dynamics can be used to derive the canonical theory of the incompressible fluid following the Dirac 
prescription. The application of this theory will be discussed in following publication. 
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Appendix A 

Matrix M: 

The matrix elements M ab satisfy the following system of partial differential equations: 



V 7 



V z M 2 i(x,z) + _Lv z M 22 (x,z) x (V x v(z)) = 5(x-z), 
P( z ) 

V x Mi 2 (x, z) + - TT VxM 22 (x, z) x (V x v(x)) = 5(x - z), 



(15) 
(16) 
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A z Mi 2 (x,z) = -S(x-z), 
A x M 21 (x,z) = <5(x-z), 



V z • V z Afn(x,z) + — -V z Mi 2 (x,z) x (V x v(z)) 
P( z ) 



V x - V x M n (x,z) + — -V x M 21 (x,z) x (V x v(x)) =0, 
A x ) 

A z Af 22 (x,z) =0, 
A x Af 22 (x,z) =0. 

It is easy to check that these equations are satisfied by matrix elements given below 

1 



-V x -G(x' - z) x (V x v(x')) 



= M{G} 



M u (x,z) = - / dx'G(x - x')V x , 

M 12 (x, z) = -M 2 i(x, z) = -G(x - z) , 
M 22 (x,z) = 0, 

Details of the Dirac brackets evaluation for the ideal fluid: 

Consider the Hamiltonian (JTJ) , the Dirac bracket C p(~x),H □ reads: 

Cp(x),tf □ = {p(x),ff}-^ y"d Zl dz 2 {p(x),8 l (z 1 )}A^(z 1 ,z 2 ){e,(z 2 ), J ff} 



= V x • J(x) - J dz, A x A/ 2 i(x, z) ^V z • J(z) 

+ A x Af 22 (x, z) [V z • (v(z) x (V x v(z))) - A z (/x(v, g))] = 



(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



(23) 



(24) 



One sees immediately that the right hand side of (|24|) vanishes due to (|18|22|23[) . Here J = gv 
denotes the fluid particle current and p(v,g) = |v| 2 /2 + df(p(z))/dg is the moving fluid chemical 
potential. 

The continuity equation is then 



d_ 

Of 



p(x,i) =C p(x,t),H^= 



(25) 



as expected. 

Evaluating Dirac bracket □ u*(x), H □ we obtain: 



C^(x),ffD = {v l ( X ),H}-J2 J {u i (x),e o (z 1 )}M o6 (z 1)Z2 ){e 6 ( Z2 ),JJ} 

a, 6 

= A^-^A^- (26) 
After straightforward but lengthy calculations we obtain: 
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A 1 



K(x),ff} = [v(x) x (V x v(x))r - V 1 [ M (v,e)] , 
-V* /dzM u (x,z) [V z - J(z)l , 



^22 = 



.4 



12 



Vj, /" dz G(x - z) {V z • [v(z) x (V x v(z))] - A z [ M (v, g)]} , 



A l 21 = - /cfe — [V x G(x-z) x (V x v(x))] 1 V z -J(z 
p(x) L 



Using above, together with equations (|T7|) and (j2~Tj) we obtain: 



-v(x,t) = v(x) x (V x v(x)) 



(27) 



V x Mn(x, z) + — -V x G(x - z) x (V x v(x)) 
/" dz [V x G(x - z)] {V z • [v(z) x (V x v(z))]} . 



V z • J(z) 



(28) 



Acting on both sides of equation ([28} with operator div = V x -, and using equations (|16|20j) one 



gets: 



- [V x • v(x, t)} =C V x • v(x, i), fT □= V x - E v(x, t), □= V x ■ -^(x, t) = 0. (29) 



Thus O2 (x) = V x ■ v(x, t) is a constant of motions, as expected. Condition p — po implies that 
V x ■ J(x, t) is also a constant of motions, and V x • </(x, t) = 0. 

Now, using equations (|23I29I) one easily sees that equation (f28|) reduces to 



9v(x,t) 
9t 



v(x,t) x (V x v(x,i)) - V x 



y dz G(x - z) V z • [v(z, *) x (V x v(z, *))] 



(30) 



which is exactly the Euler equation for an ideal, incompressible fluid in its integral form. 
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